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Abstract. We study theoretically interaction potentials and low energy collisions between 
different alkali atoms and alkali ions. Specifically, we consider systems like X + Y+, where 
X(Y+) is either Li(Cs+) or Cs(Li+), Na(Cs+) or Cs(Na+) and Li(Rb+) or Rb(Li+). We 
calculate the molecular potentials of the ground and first two excited states of these three 
systems using pseudopotential method and compare our results with those obtained by others. 
We derive ground-state scattering wave functions and analyze cold collisional properties of 
these systems for a wide range of energies. We find that, in order to get convergent results for 
the total scattering cross sections for energies of the order 1 K, one needs to take into account 
at least 60 partial waves. Low energy scattering properties calculated in this paper may serve 
as a precursor for experimental exploration of quantum collisions between an alkali atom and 
an alkali ion of different nucleus. 
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l. Introduction 

Over more than last three deeades, there has been tremendous progress in the teehnology of 
eooling and trapping of single-eleetron alkali atomie gases and some two-eleetron atomie 
gases. In parallel, the developments in the teehnology of trapping and cooling of some 
ions, specially alkaline-earth cations over last several decades has led to new vistas of 
research activities with trapped and cold ions. With recent experimental developments 
|[Ill2l[3l|4l[51[6l|71[8l|9l [lol [IH [I2l [T3l [Ml [T5l [l6|| using hybrid ion-atom traps [[I71|, both 
atoms and ions can be confined simultaneously in a common space or an atom-trap can be 
merged with an ion-trap. These hybrid systems open up prospects for emulating solid-state 
physics with laser-generated periodic structure of ions ifT^ and for exploring physics of ion- 
controlled Josephson junction |[T9l l20]| . Tonk-Girardeau gas with an ionic density “bubble” 

m. 

At a fundamental level, these new hybrid devices facilitate the laboratory investigations 
on ion-atom collisions in hitherto unexplored low energy regimes - from milliKelvin 
down to microKelvin or sub-microKelvin temperature regime. Ion-atom scattering at low 
energy is important for a number of physical systems, such as cold plasma, planetary 
atmospheres, interstellar clouds. Gaining insight into the ion-atom interactions and scattering 
at ultra-low energy down to Wigner threshold law regime is important to understand charge 
transport [l22l at low temperature, radiative association |l2l[3l|23l, ion-atom bound states (Mi, 
ion-atom photoassociation ll^ [26i . and many other related phenomena. Several theoretical 
investigations (271 l28l l29l [30l (SB [3B EB |34l of atom-ion cold collisions have been carried 
out in recent times. It is proposed that controlled ion-atom cold collision can be utilized 
for quantum information processes (35l [36i. Several recent theoretical studies have focused 
on charge-transfer (371 EH [39i and chemical reaction processes (401 mi at ultralow energies 
between an alkaline-earth ion or Yb+ and an alkali atom. 

Ions are usually trapped by radio-frequency fields. The major hindrance to cool trapped 
ions below milliKelvin temperature stems from the trap-induced micro-motion of the ions, 
which seems to be indispensable for such ion traps. Therefore, it is difficult to achieve 
sub-microKelvin temperature for an ion-atom system in a hybrid trap. But, to explore fully 
quantum or Wigner threshold law regime for ion-atom collisions, it is essential to reduce 
the temperature of ion-atom hybrid systems below one micro-Kelvin. One way to overcome 
this difficulty is to device new methods to trap and cool both ions and atoms in an optical 
trap. With recent experimental demonstration of an optical trap for ions (42l . the prospect for 
experimental explorations of ultracold ion-atom collisions in Wigner threshold regime appears 
to be promising. 

In the contexts of ion-atom cold collisions of current interest, there are mainly four types 
of ion-atom systems:(l) an alkaline-earth ion interacting with an alkali metal atom, (2) a rare 
earth ion (such as Yb) interacting with an alkali metal atom, (3) an alkali metal ion interacting 
with an alkali atom of the same nucleus and (4) an alkali metal ion interacting with an alkali 
atom of different nucleus. In the first two types, both the ion and the atom have one valence 
electron in the outermost shell. In the third type, the ion has closed shell structure while the 
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Figure 1. Adiabatic potential curves for (LiRb)+, (LiCs)+ and (NaCs)+ 


atom has one valence electron in the outermost shell. Since both the ion and atom are of the 
same nucleus, there is center of symmetry for the electronic wave function of the ion-atom 
pair. Because of this symmetry, resonant charge transfer collision is possible in the third type. 
In the fourth type, the ion is of closed shell structure while the atom has one valence electron 
in the outermost shell. Since the nuclei of the atom and the ion are different, there is no res¬ 
onant charge transfer collision in fourth type, though non-resonant charge transfer is possible 
in all the types. For scattering in the ground-state potentials, non-resonant charge transfer 
collision is suppressed in the ultralow energy regime (< /xK). 

Here we study cold collision in the fourth type of ion-atom systems, though these 
systems are not straightforward to obtain experimentally. Normally, neutral alkali-metal 
atoms and alkaline-earth ions can be laser-cooled and trapped. However, there are some recent 
experimental studies showing different techniques ifTOlfTTlfldlfTSll for trapping and cooling of 
closed-shell akali-metal ions. Our primary aim is to understand ground-state elastic scattering 
processes between such an ion and an alkali atom of different nuclei at ultralow energy where 
inelastic charge transfer process can be ignored. A proper understanding of ground-state 
elastic scattering processes (SSI [34l| at low energy in an ion-atom system is important for 
exploring coherent control of ion-atom systems and formation of ultracold molecular ions by 
radiative processes. Due to the absence of resonance charge transfer reaction, the fourth type 
system is preferable for our purpose. However, there remains finite possibility for inelastic 
collisions if the hyperfine structure of the neutral atom is taken in account. In our calculations 
we do not consider hyperfine interactions. 

Let a cold alkali metal ion A+ interact with a cold alkali metal atom B. Suppose, the 
atomic mass of B is smaller than that of atom A. Then, for AB+ molecular system, ground- 
state continuum asymptotically corresponds to the atom B (ns) in electronic state and the 
ion A+ (complete shell) in electronic state. In the separated atom limit, the first excited 
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molecular potential asymptotically goes as the atom A (ns) in electronic state and alkali 
ion (complete shell) in electronic state while the second electronic excited potential 
corresponds to B (np) in electronic state and alkali ion A"*" (complete shell) in electronic 
state. 


Table 1. A comparison of the spectroscopic constants for the ground (X^S+) and the first and 
second excited (2^E+ and 3^S+) electronic states of (LiRb)+ molecular ion with the available 
works. 


State 

i?e(a.u.) 

L)e(cm ^) 

Te(cm ^) 

uieicm ^) 

uJeXeicm 1) 

5e(cm ^) 

references 

X2S+ 

7.70 

3912 


137.56 

1.67 

0.158098 

This work 



3705 





ttH 


7.50 






[El 


7.60 

3432 





lEl 


7.54 

4193 


139.65 



lEa 

22S+ 

12.67 

1270 

12438 

63.73 

0.51 

0.058456 

this work 


12.60 

1306 


63.00 



lEa 

32S+ 

16.66 

617 

18200 

35.73 

1.55 

0.033769 

This work 


We investigate cold collisions over a wide range of energies in three different alkali 
atom-alkali ion systems in the ground molecular potentials. The colliding atom-ion pairs we 
consider are: (i) Cs’*' -i- Li, (ii) Cs"*" -i- Na, and (iii) Rb+ -i- Li. Since at ultralow energy, charge 
transfer reaction is highly suppressed in these systems, we study elastic scattering only. To 
calculate scattering wave functions, we need the data for Born-Oppenheimer adiabatic poten¬ 
tials of the systems. In the long range where the separation r > 20 oq (oq is the Bohr radius), 
the potential is given by the sum of the dispersion terms, which in the leading order goes 
as 1/r^. We obtain short-range potentials by pseudopotential method. The short-range part 
is smoothly combined with the long-range part to obtain the potential for the entire range. 
We then solve time-independent Schroedinger equation for these potentials with scattering 
boundary conditions by Numerov-Cooley lf4^ algorithm. We present detailed results of scat¬ 
tering cross sections for these three systems for energies ranging from 0.01 micro-Kelvin 
(/iK) to 1 Kelvin (K). However, the low energy regime of our interest ranges from sub-/iK to 
1 milli-Kelvin (mK). Here we choose to use Kelvin as the unit of energy to readily convey 
the information how cold the system should be to explore such collision physics. Normally 
the results of atomic collisions are expressed in atomic unit (a.u.). Note that 1 a.u. of energy 
corresponds to 3.1577465 x 10^ K or 27.21138386 electron-Volt (ev). 
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Figure 2. (Color online) In the upper panel, the centrifugal energies in unit of milliKelvin 
(mK) for s- (black, solid), p- (red, dashed) and d- (green, dashed-dotted) partial waves are 
plotted against atom-ion separation for 1^1]+ state of (LiCs)+. The lower panel shows the 
same for (NaCs)+. 


2. Interactions and molecular properties 

Under Born-Oppenheimer approximation, we eompute the adiabatie potential energy eurves 
of the and eleetronie states of the three ionic molecules (LiRb)’*', (LiCs)’*' 

and (NaCs)+ using the pseudopotential method proposed by Barthelat and Durand ll4^ in 
its semi-local form and used in many previous works [|431 14^ l47l l48l 1^ . The interaction 
potentials between different alkali metal ion and an alkali metal atom were calculated by 
Valance ll50]l in 1978 by the method of pseudopotential. Here we briefly describe spectro¬ 
scopic constants and interaction potentials of the three considered atom-ion systems, and the 
pseudopotential method used to obtain these molecular properties. Our final goal is to calcu¬ 
late cold collisional properties of these systems using these potentials as described in the next 
section. 


2.1. Adiabatic Potentials 

For the three ionic systems (LiRb)+, (LiCs)+ and (NaCs)+, the potential energy curves of the 
1-3^S+ electronic states, are built from an ab initio calculation for the internuclear distances 
ranging from 3 to 200 Oq (oq is Bohr radius) as discussed in the sub-section 2.2. These states 
dissociate respectively, into Li(2s and 2p) -i- Rb+ and Rb(5s) -i- Li+, Li(2s and 2p) -i- Cs’*' 
and Cs(6s) -i- Li+, and Na(3s and 3p) -i- Cs’*' and Cs(6s) -i- Na+. They are displayed in figure 
1(a) for LiRb+, 1(b) LiCs’*' and 1(c) for NaCs"*". In the long range part beyond 200ao, the 
potential is given by the sum of the dispersion terms. The short range part (< 20ao) of ab 
initio potential is smoothly combined with the long-range part by interpolation to obtain the 
potential for entire range extending to several thousand oq. For the three ionic systems, we 
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Figure 3. (Color online) s- p- and d-wave scattering cross sections in atomic unit (a.u.) or in 
unit of Oq (ao is Bohr radius) are plotted against collision energy E in Kelvin (K) for 
state of (LiCs)+ and (NaCs)+ in upper and lower panel, respectively. 


remark, that the ground state has the deepest well compared to the and excited 
states. Their dissociation energies are of the order of several 1000 cm~^. Their equilibrium 
positions lie at separations that are relatively larger than those of typical neutral alkali-alkali 
diatomic molecules. 


The long range potential is given by the expression 


V{r) 


1 

2 



( 1 ) 


where C^, Cq correspond to dipole, quadrupole polarisabilities of concerned atom. Hence, 
the long range interaction is predominately governed by polarisation interaction. Dipole 
polarisabilities for Na (3s) and Li (2s) are 162 a.u. and 164.14 a.u., respectively. One can 
define a characteristic length scale of the long-range potentials hy j3 = The 

values of fi for the collision of Li - Cs^ pair and Na - Cs'*' pair are 1411.5 a.u. and 2405.8 
a.u., respectively 

For the ground state l^S'*' the well depth and the equilibrium position are 2979 cm“^ 
and 8.51 a.u., respectively, for (NaCs)'*’ while those for (LiCs)'*' are 3176 cm“^ and 8.12 a.u., 
respectively. Reduced mass for (NaCs)"*" and (LiCs)'*' (for ^Li isotope) are taken as 19.5995 
a.u. and 6.6642 a.u., respectively. For both ®Li-Rb+ and ®Li-Rb+ isotopes, we used the 
reduced masses 5.6171 and 6.4805 a.u. respectively. 


2.2. Pseudopotential method: Results and discussions 

The use of the pseudopotential method, for each (XY)+ ionic molecules (LiRb+, LiCs’*' and 
NaCs'*'), reduce the number of active electrons to only one electron. We have used a core 
polarization potential Vcpp for the simulation of the interaction between the polarizable X"*" 
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Figure 4. (Color online) Logarithm of total scattering cross section in a.u. as a function of 
logarithm of energy in K for electronic ground-state 1^E+ of (LiCs)+ and (NaCs)+ are 
plotted in upper and lower pannel, respectively. 


and Y+ cores with the valenee eleetron. This eore polarization potential is used aeeording to 
the formulation of Muller et al. lBTll . and is given by 

Vcpp = —-SaGa/a • fx (2) 

where ctA and fx are respeetively the dipole polarizability of the eore A, and the eleetrie field 
produced by valenee eleetrons and all other eores on the core A. The eleetrie field fx is defined 
as: 

fx = ^^'^F{fa,px) - ^y^x^Zx (3) 

^y,\ 

where Vix and Ry,x are respeetively the eore-eleetron veetor and the core-core veetor. 

Based on the formulation of Fouerault et al. l[52]| the eut-off funetion F(rjA, Pa) is a fune- 
tion of the quantum number 1. In this formulation, the interaetions of valenee electrons of 
different spatial symmetry with eore electrons are considered in a different way. The used 
eut-off radii of the lowest valenee s, p, d and f one-eleetron for the Li, Na, Rb and Cs 
atoms are taken from ref. [l43l |471 Il8l |49l. The extended Gaussian-type basis sets used for 
Li, Na, Rb and Cs atoms are respeetively (9s8p5dlf/8s6p3dlf) ll48l . (7s6p5d3f/6s5p4d2f) ll48ll . 
(7s4p5d 1 f/6s4p4d 1 f) ifS^ and (7s4p5d 1 f/6s4p4d 1 f) Hbll . The used eore polarizabilities of Li+, 
Na+, Rb+ and Cs+ eores, whieh are equal respeetively 0.1915 Oq, 0.993 Oq, 9.245 Oq and 
15.117 Oq are taken from llSTI . Using the pseudopotential teehnique eaeh moleeule is redueed 
to only one valenee eleetron interaeting with two eores. Within the Born-Oppenheimer ap¬ 
proximation an SCF calculation, provide us with accurate potential energy curves and dipole 
functions. 
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2.3. Spectroscopic constants 


The spectroscopic constants {Re, De, Tg, ccg, oJeXe, Be) of the 1-3^S+ electronic states are 
presented in tables 1-3 for (LiRb)+, (LiCs)’*' and (NaCs)’*' respectively. To the best of our 
knowledge, no experimental data has been found for these systems. We compare our spectro¬ 
scopic constants only with the available theoretical results. Table 1 presents the spectroscopic 
constants of (LiRb)+ compared with the other theoretical works [l54l l55l 15^ ISTll . As it seems 
from table 1, there is a good agreement between the available theoretical works [l54l[55l 15^ISTll 
and our ab initio study. Azizi et al. llSTll reported the spectroscopic constants for many ionic 
alkali dimers but only for the ground and the first excited states. They have used in their 
study a similar formalism as used in our work. We remark that our ground state equilibrium 
distance {Re) presents a satisfying agreement as well as the well depth {De) with the work of 
Azizi et al. llSTll . We found for {Re) and {De), respectively, 7.70 a. u. and 3912 cm“^ and they 
found 7.54 a.u. and 4193 cm“^. The difference between the results of Azizi et al. [57] and 
our values are about 2.12% and 6.70% for Re and De, respectively. The same agreement is 
observed between our harmonacity frequency (ccg) and that of Azizi et al. lfSTll . The difference 
between the two values is 2.09 cm“\ which represents a difference of 1.49%. There is also a 
very good agreement between our equilibrium distance and that of Patil et al. [[5^ {Re= 7.60 
a.u.), in contrast to their well depth lfS^ . which is underestimated {De=34-32 cm~^). For the 
first excited state, Azizi et al. lfSTll reported the spectroscopic constants i?e=12.60 a.u., De = 
1306 cm“^ and Ue = 63.00 cm~^ to be compared with our values of, respectively, 12.67 a.u., 
1270 and 63.73 cm”k The differences in percentages between the values obtained by Azizi 
et al. [57] and our results are 0.55%, 2.75% and 1.14% for Re, De and Ue, respectively. 


The spectroscopic constants of (LiCs)+ and (NaCs)+, are presented respectively in tables 
2 and 3 and compared with the available theoretical results ll58l [50]l . To the best of our 
knowledge there is no experimental data for these ionic molecules. For the (LiCs)+ ionic 
molecule, we compare our results only with the theoretical work of Khelifi et al. ll5^ where 
they have used a similar method as used in our work. They reported for the ground state 1^S+ 
the spectroscopic constants Re = 7.91 a.u. and De = 3543 cm~^ to be compared with our 
values of, respectively. Re = 8.12 a.u. and De = 3176 cm“^. A rather good agreement is 
observed for the equilibrium distance, however their potential is much deeper. Reasonably 
good agreement between our spectroscopic constant and those of Khelifi et al. [|5^ is 
observed for the 2^S+ and 3^S+ states. The difference between Khelifi et al. spectroscopic 
constants and our results are found to be 2.58% and 10.35% for Re and De, respectively. For 
the ground state of (NaCs)’*' ionic molecule, there is a good agreement between our well- 
depth, as well as the equilibrium distance with the theoretical results of Valance [|50]| . We 
found a well-depth of 2979 cm"^ located at 8.51 a u., while valance ll50ll found 3388 cm~^ 
located at 7.60 a.u.. The difference in percentages between valances data and our values are 
12.07% and 11.97% for Re and De, respectively. In contrast to the ground state, where the 
agreement between our work and those of Valance lISOll is good, there is a disagreement for the 
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Figure 5. (Color online) Same as in figure 4 but for Li + (1^E+) collision. The black 

solid line and red dotted lines are for ®Li and ^Li, respectively. Black dashed lines represent 
linear fit to both the curves for energies greater than 10“® K. Both cases follow one-third law 
and the intercepts are found to be 3.587 and 3.606 in a.u. for ®Li and ^Li, respectively; where 
as the theoretically calculated values are 3.3269 and 3.3476 in a.u., respectively. 


and states. The two states exhibit small potential wells in Valances work equal, 
respectively, to 774 and 732 cm~^; while in our work we have found significantly higher well 
depths of 1262 and 1393 cm“^, respectively. Our well depth for the 2^S+ and 3^E+ states are 
of 63.04% and 90.30% larger than those of Valance [50]. 

The large discrepancies between the Valences results and our values can be explained by 
two reason. First, Valence used the Hellmann-type model where the Hamiltonian is written 
in terms of one valence electron interacting with closed shell cores. Therefore, the core-core 
interaction potentials were neglected. This is certainly true for large interatomic distances 
but it is not accurate enough for small separations. Valence considered this shortcoming 
as not relevant to the charge exchange collision process they studied using the interaction 
potential where only intermediate and long-range distances are required. Although, without 
considering the core-core interaction one can get correct asymptotic limits, at intermediate 
distance this will affect the accuracy of the equilibrium distance and well depth. Second, 
the core-valence interactions were considered in our work using a formulation where, for 
each atom, the core polarization effects are described by an effective potential as described 
previously. The latter is modified by an I- dependent cut-off function F to consider, in a 
different way the orbital symmetries. In the case of Valence, the Hellmann pseudopotential is 
not dependent and the parameters are optimized to obtain the two first experimental energy 
levels. Similarly, in the the absence of the core-core interaction. Valence considered that this 
lack of ’7” dependence is not crucial for the charge exchange collision. In our approach, 
in addition to the use of an /-dependent pseudoptential the parameters were optimized to 
reproduce with high accuracy the ionization potential and the experimental energy levels of 
many alkali excited states. Moreover, in our study, the interatomic separation grid is more 
extended as R varies from 3 to 200 a.u. with a distance step of 0.1 a.u, but Valence limited 





Interactions and low energy collisions between an alkali ion and an alkali atom of different nucleuslO 

his calculation for intemuclear distances ranging from 5 a.u. to 25 a.u. It is important 
to note that the ealeulation time in our ealeulation is enormously redueed due to the use 
of the pseudopotential approaeh that leads these ionie moleeules to one-eleetron systems. 
Caleulations are performed on a desktop Pentium IV Aeer eomputer where the time for one 
single distanee is about 32 seeonds. 

3. Elastic collisions: Results and discussions 

Here we present results on elastie eollisions of the diseussed ion-atom systems by FORTRAN 
eode [l59ll using the well-known Numerov-Cooley algorithm [|43]|. This algorithm provides an 
effieient teehnique for numerieally solving a seeond order differential equation. It uses a three- 
point reeursion relation to ealeulate first order derivative and makes use of exaet seeond order 
derivative provided by the equation itself. To solve time-independent Sehroedinger equation 
as a seattering problem, one has to propagate the Numerov-Cooley eode from initially small 
intemuelear separation r to asymptotieally large r where the wave funetion eorresponds to the 
state of a free partiele and so behaves sinusoidally. The seattering S'-matrix is then dedueed 
by matehing the asymptotie solution with the standard seattering boundary eonditions. The 
eorreetness of the ealeulations is ensured by the unitarity of the S'-matrix. The lower the 
energy larger is the asymptotie boundary. To initiate the propagation of the eode from r ~ 0 
position, one has to set initially the values of the wave funetions at two initial positions, then 
the eode ealeulates the wave funetion for the third position, and then taking the the values 
of the wave funetions for the seeond and third position it ealeulates the wave funetion for the 
fourth position, and the proeess eontinues until the asymptotie boundary is reaehed. The initial 
boundary eonditions are set by expanding the interaetion potential for small r and solving the 
Sehroedinger equation analytieally for r —)■ 0 limit. A good measure for the asymptotie 
seattering boundary (large r) ean be found by setting the eondition that the effeetive potential 
Veff = —Cq/t^ + (/i^/2/i)£(£ -f l)/r^ for large r beeomes mueh less than the eollision energy 
E (at least one tenth of E). 

We use the interaetion potentials ealeulated above as input data. In partieular, we focus 
on low energy eollisions. Using well-known expansion of eontinuum state in terms of partial 
waves, the wave funetion ( k r) for £th partial wave is given by 



(4) 


subjeet to the standard seattering boundary eondition 

'fe,{kr) ~ sin [kr — in/2 + rji] 


(5) 


Here r denotes the ion-atom separation, the wave number k is related to the eollision energy 


Ehy E = b? fc^/2/i and /r stands for the redueed mass of the ion-atom pair. The total elastie 
seattering eross seetion is given by 



CO 
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To know the relevant energy regimes where s- p- and d-wave collisions are important, 
we have plotted in figure 2 the centrifugal energies for first 3 partial waves (£ = 0, 1 and 2) 
against ion-atom separations for ground state 1^S+ of (NaCs)’*' and (LiCs)’*', respectively. 

For d-wave, the values of centrifugal barrier are about 0.007 mK for (NaCs)+ and about 0.06 
mK for (LiCs)+. These values indicate that the potential energy barriers for low lying higher 
partial waves are very low for atom-ion systems allowing tunneling of the wave function 
towards the inner region of the barriers. Unlike atom-atom systems at low energy, a number 
of partial waves can significantly contribute to the ion-atom scattering cross section at low 
energy. We find that in order to get convergent results at milli- and micro-Kelvin regimes, the 
numerical calculation of scattering wave function needs to be extended to at least 10000 oq 
and 20000ao, respectively. . 

In figure 3, the partial-wave cross sections for ground 1^S+ state collisions are plotted 
against collision energy U in K for (NaCs)’*' and (LiCs)"'". As A; —)■ 0 , s-wave cross section 
becomes independent of energy while p- and d-wave cross sections vary in accordance with 
Wigner threshold laws. According to Wigner threshold laws, as k ^ 0, the phase shift for 
Ah partial-wave pi ~ if £ < (n — 3)/2, otherwise pi ~ for a long range potential 
behaving as l/r"^. Since Born-Oppenheimer ion-atom potential goes as 1/r^ in the asymptotic 
limit, as fc —> 0, s-wave scattering cross section becomes independent of energy while all 
other higher partial-wave cross sections go as As energy increases beyond the Wigner 
threshold law regime, the s-wave cross section exhibits a minimum at a low energy which 
may be related to Ramsauer-Townsend effect. It can be noticed from figure 3 that at Ramsauer 
minimum of s-wave scattering, p- and d-wawe elastic cross sections are finite. This means that 
one can explore p- or d-wave interactions of ion-atom systems at the minimum position that 
occurs at a relatively low energy (/i- and milli-Kelvin regime). 

In the Wigner threshold regime, hetero-nuclear ion-atom collisions are dominated by 
elastic scattering processes since the charge transfer reactions are highly suppressed in this 
energy regime. Furthermore, resonant charge transfer collisions do not arise in collision of an 
ion with an atom of different nucleus. 

It is worthwhile to compare the elastic scattering results of hetero-nuclear alkai ion-atom 
system to those of homo-nuclear alkali atom-ion system Na-i-Na+ ll22]l . By comparing our 
results in figure 3 with the figures 2 and 3 of Ref. Il22ll . we notice that elastic scattering 
cross sections of both homo- and hetero-nuclear ion-atom systems at ultralow energies are 
comparable and primarily governed by Winger threshold laws. However, unlike hetero- 
nuclear systems, the resonant charge transfer cross section in homo-nuclear Na-i-Na+ is also 
of the same order as that of elastic one at ultra-low energies. Resonant charge exchange cold 
collisions between Yb and Yb+ are experimentally observed and found to be a dominant 
process in this homo-nuclear system [[U. 

In figure 4, we have plotted logarithm of total elastic scattering cross sections (atot) 
expressed in a.u. as a function of logarithm of U in K for ground state collisions between Li 
and Cs+ and between Na and Cs+, respectively. For both (LiCs)’*' and (NaCs)+, to calculate 
a tot, we require more than 60 partial waves to get converging results for energies greater than 
1 mK. This is because, as the collision energy increases, more and more partial waves start to 


nd an alkali atom of different nucleus\2 



Figure 6. (Color online)Partial wave cross sections for ®Li + Rb+ (1^I]+) and ^Li + Rb+ 
(1^I]+) collision are plotted as a function of E (in K) for i = 0 (solid) and £= I (dashed) in 
upper and lower panel, respectively. 


contribute to a tot- In order for a partieular partial wave i to eontribute appreeiably at a given 
energy E, the centrifugal barrier height eorresponding to this partial wave should not be too 
high eompared to the energy. An estimate of the minimum number of partial waves required 
for obtaining eonvergent results may be made by taking the height of the eentrifugal barrier to 
be higher than the energy by one order of magnitude. In this context, it is worth mentioning 
that the eentrifugal barrier height for a partial wave of an ion-atom system is usually mueh 
lower than the eorresponding height of atom-atom system, therefore larger number of partial 
waves contribute to atot in ease of an ion-atom system. In the large energy limit, the cross 
section behaves as 


(^tot ~ TT 




(7) 


Thus, in the large energy regime the slope of the logarithm of atot as a function of log E 
is a straight line obeying the equation atot{E) = —(l/3)i? -|- ce where the slope of the line is 
— 1/3 and the intereept of the line on the energy axis is solely determined by the C 4 eoeffieient 
of long-range part of the potential, or equivalently by the characteristic length scale (3 of the 
potentials or the polarizability of the neutral atom interaeting with the ion. As shown in 
figure 4, we have numerieally verified this law for both Li-Cs^ and Na-Cs'*' systems. 
The linear fit to the the plot of numerically calculated atot against in the large energy regime 
shows that the value of the slope is quite close to the aetual value of 1/3. Sinee the dipole 
polarizability of Li and Na is not much different, one can expect that for both Li-Cs"*" and 
Na-Cs"*" systems, the energy-dependenee of atot at large energy should be similar. In fact, 
figure 4 shows that for energies much greater I /iK, both systems exhibit similar asymptotie 
energy dependenee. 

Finally, we consider the eollision of ®^Rb+ with the two isotopes of Li atom, namely, ®Li 
and ^Li, and the results are shown in figures 5 and 6 . The purpose here is to investigate the 
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isotopic effects of Li on low energy collisions with ion . In figure 5 we have plotted 

logarithm (to the base of 10) of a tot as a funetion of logarithm (to the base of 10) of E. Results 
elearly show that the patterns are same for ®Li and ^Li eolliding with ®^Rb+, but they differ in 
magnitudes, partieularly in the low energy regime. Figure 6 exhibits the variation of s- and p- 
wave seattering eross seetions against energy in log-log seale for both isotopes of Li eolliding 
with Rb+. Again, the results are qualitatively similar for both the isotopes, but they differ 
slightly quantitatively. Comparing the upper and lower panels in figure 6, one ean notice that 
the Wigner threshold law regime for ®Li can be attained at slightly higher energy than that for 
"Li. 

4. Conclusions 

In eonelusion, we have studied elastie seattering between alkali ion and alkali atom of different 
nuclei of three ion-atom systems, namely, Li -i- Cs"*", Na -i- Cs+, and Li -i- Rb+. We have 
ealeulated the interaetion potentials and spectroseopie eonstants of these systems. We have 
presented a detailed study of elastic collision physics over a wide range of energies, showing 
the onset of Wigner threshold regime at ultra low energy and the 1/3 law of seattering at higher 
energy regime. The low energy seattering results presented here may be useful for future 
exploration for radiative- or photo-assoeiative formation of eold diatomic molecular ions from 
these three ion-atom systems. The colliding ground state atom-ion pair corresponding to the 
eontinuum of l^E^ may be photoassoeiated to form a bound state in second excited potential 
3^S+ by using a laser of appropriate frequeney. This is promising beeause the transition is 
dipole allowed and is similar to atomie transition. The formation of sueh eold moleeular ions 
by radiative proeesses will pave new directions to ultracold chemistry. 
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Table 2. A comparison of the spectroscopic constants for the ground (X^S+) and the first and 
second excited (2^S+ and 3^S+) electronic states of (LiCs)+ molecular ion with the work of 
Khelifi et al.^^ 


State 

i?e(a.u.) 

£>e(cm 1 ) 

Te(cm 1) 

uie(cm 

oJeXeicm 1) 

Ee(cm 

references 

12 S+ 

8.12 

3176 

0 

124.46 

1.08 

0.138337 

This work 


7.19 

3543 





[ISl 

22 S+ 

12.42 

1911 

13343 

68.69 

0.63 

0.059126 

This work 


12.37 

2022 

19553 




[ISl 

32S+ 

18.53 

422 

17660 

28.07 

1.29 

0.026570 

This work 


18.38 

409 

20849 




IlSl 


Table 3. A comparison of the spectroscopic constants for the ground and the first and 

second excited (2^E+ and electronic states of (NaCs)'*' molecular ion with the work of 

Valance'^^. 


State 

i?e(a.u.) 

De{cm 1) 

Te(cm 1) 

ujeiom 1) 

OJeXeicm 

5e(cm 

references 

12E+ 

8.51 

2979 

0 

68.17 

0.32 

0.042418 

This work 


7.60 

3388 





IlSoli 

22S+ 

13.45 

1262 

11758 

33.98 

0.30 

0.016977 

This work 


14.20 

774 





IlSoli 

32S+ 

16.77 

1393 

18553 

26.20 

1.69 

0.010920 

This work 


15.00 

732 





ttsa 








